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We present a complete analysis of multipartite entanglement of three-mode Gaussian states of continuous 
variable systems. We derive Standard forms which characterize the covariance matrix of pure and mixed three- 
mode Gaussian states up to local unitary operations, showing that the local entropies of pure Gaussian states are 
bound to fulfill a relationship which is stricter than the general Araki-Lieb inequality. Quantum correlations can 
be quantified by a proper convex roof extension of the squared logarithmic negativity, the continuous-variable 
tangle, or contangle. We review and elucidate in detail the proof that in multimode Gaussian states the contangle 
satisfies a monogamy inequality constraint [G. Adesso and F. Illuminati, New J. Phys. 8, 15 (2006)]. The 
residual contangle, emerging from the monogamy inequality, is an entanglement monotone under Gaussian local 
operations and classical communication and defines a measure of genuine tripartite entanglement. We determine 
the analytical expression of the residual contangle for arbitrary pure three-mode Gaussian states and study in 
detail the distribution of quantum correlations in such states. This analysis yields that pure, symmetric states 
allow for a promiscuous entanglement sharing, having both maximum tripartite entanglement and maximum 
couplewise entanglement between any pair of modes. We thus name these states GHZ/VF states of continuous 
variable systems because they are simultaneous continuous-variable counterparts of both the GHZ and the W 
states of three qubits. We finally consider the effect of decoherence on three-mode Gaussian states, studying 
the decay of the residual contangle. The GHZ/W states are shown to be maximally robust against losses and 
thermal noise. 

PACS numbers: 03.67 .Mn, 03.65.Ud 



I. INTRODUCTION 

Multipartite entanglement is one of the most fundamental 
and puzzling aspects of quantum mechanics and quantum in- 
formation theory. Although some progress has been recently 
gained in the understanding of the subject, many bàsic prob- 
lems are left to investigate in this fascinating area of research. 
Multipartite entanglement poses a bàsic challenge both for the 
obvious reason that it is ubiquitous to any practical realization 
of quantum communication protocols and quantum computa- 
tion algorithms, and because of its inherent, far-reaching fun- 
damental interest 001. 

The steps undertaken so far in the attempt to reach some un- 
derstanding of quantum entanglement in multipartite settings 
can be roughly classified in two categories. On the one hand, 
the qualitative characterization of multipartite entanglement 
can be investigated exploring the possibility of transforming a 
multipartite state into another under different classes of local 
transformations and introducing distinct equivalence classes 
of multipartite entangled states |2|. On the other hand, a quan- 
titative characterization of the entanglement of states shared 
by many parties can be attempted: this approach has lead to 
the discovery of so-called monogamy inequalities, constrain- 
ing the maximal entanglement shared by different internal par- 
titions of a multipartite state |3|, |4{] . Such inequalities are up- 
rising as one of the possible fundamental guidelines on which 
proper measures of multipartite entanglement should be built. 



Recently, much effort has been devoted to the study of en- 
tanglement in continuous-variable systems, focusing both on 
quantum communicationprotocols and on fundamental the- 
oretical issues @> El- A rich and complex structure 
has emerged, already in the restricted, but physically relevant, 
context of Gaussian states. The genèric study of Gaussian 
states presents many interesting and appealing features, be- 
cause it can be carried out exploiting the powerful formalism 
based on covariance matrices and symplectic analysis. These 
properties allow to face and answer qüestions that are in gen- 
eral much harder to discuss in discrete variable systems, and 
open up the possibility to shed some light upon general facets 
of multipartite entanglement, that might carry over to systems 
of qubits and qudits. 

For two-mode Gaussian states, the qualification and quan- 
tification of bipartite entanglement have been intensively stud- 
ied, and a rather complete and coherent understanding begins 
to emerge l^. fïoll . However, in the case of three-mode Gaus- 
sian states, the simplest non-trivial instance of multi-party en- 
tangled Gaussian states that can be conceived, the multipar- 
tite sharing structure of quantum correlations presents several 
subtle structural aspects that need to be elucidated. There- 
fore, three-mode Gaussian states constitute an elementary but 
very useful theoretical laboratory that is needed toward the 
understanding of the patterns by which quantum correlations 
distribute themselves among many parties. 

A fairly complete qualitative characterization of entan- 
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glement in three-mode Gaussian states has been recently 
achieved Ulli . In the present paper, we study and present a 
fully quantitative characterization of entanglement in three- 
mode Gaussian states. We discuss the general properties 
of bipartite entanglement in pure and mixed states as well 
as the definition and determination of monogamy inequali- 
ties, genuine tripartite entanglement, and the ensuing struc- 
ture of entanglement sharing. We single out a special class 
of pure, symmetric, three-mode Gaussian states that are the 
continuous-variable anàlogues and possess the same entangle- 
ment properties of both the W and the Greenberger-Horne- 
Zeilinger (GHZ) maximally entangled states of three qubits. 
Finally, we discuss the decoherence of three-mode Gaussian 
states and the decay of tripartite entanglement in the presence 
of noisy environments, and outline different possible gener- 
alizations of our results to n-mode Gaussian states with arbi- 
trary n. 

The paper is organized as follows. In Section II we pro- 
vide a self-contained introduction to the symplectic formalism 
for covariance matrices, and review the structure of entangle- 
ment in two-mode Gaussian states. In Section III we apply the 
known facts on two-mode states and the symplectic formalism 
to provide a systematic quantification of bipartite entangle- 
ment in three-mode Gaussian states. In Section IV we review 
the concept of continuous-variable tangle and the continuous- 
variable monogamy inequalities recently derived and 
exploit these results to quantify the genuine tripartite entan- 
glement in three-mode Gaussian states. In Section V we an- 
alyze the distributed entanglement and the structure of en- 
tanglement sharing in three-mode Gaussian states, and iden- 
tify some classes of symmetric, pure and mixed, three-mode 
Gaussian states with special entanglement properties, includ- 
ing the so-called "GHZ/VF" states that maximize simultane- 
ously the genuine tripartite entanglement and the bipartite en- 
tanglement of any two-mode reduction. In Section VI we dis- 
cuss the decoherence of three-mode Gaussian states and the 
decay of tripartite entanglement due to the coupling with the 
environment. Finally, in Section VII we give some conclud- 
ing remarks and sketch an outlook on some future develop- 
ments and extensions to more general states and instances of 
continuous-variable systems. 



II. PRELIMINARY FACTS AND DEFINITIONS FOR 
GAUSSIAN STATES 



In this section, we will introduce bàsic facts and notation 
about Gaussian states of bosonic fields, reviewing some of 
the existing separability criteria for two-mode and multimode 
states and the computable measures of entanglement available 
for bipartite systems. Such bàsic results will be needed in 
extending the analysis to multipartite quantum correlations in 
multimode Gaussian states. 



A. Covariance matrices, symplectic eigenvalues, and 
inseparability criteria 

Let us consider a quantum system described by n pairs 
of canonically conjugated operators, for instance the quadra- 
ture operators of a bosonic field, pj}, satisfying the 
canonical commutation relations [stj, p k ] = Sjk- For ease 
of notation, let us define the vector of field operators R = 
{xi,pi, . . . , x n ,Pn} and note that the commutation relations 
can be written as [R, R] = 2ifl, where the symplectic form íl 
is defined as 

n = é". "=(_i l) > <D 

where (J) denotes the direct sum. Any state of such a system 
is represented by a hermitian, positive, trace-class operator 
g, the so-called density matrix. Gaussian states are defined 
as states with Gaussian characteristic (and quasi-probability) 
functions: a state g is Gaussian if and only if its characteristic 
function 

X (0 = Tr [qDç] , (2) 

where £ E R 2n is a real vector and Dç = cxp (iR T íl£) is 
Glauber's displacement operator, is a multivariate Gaussian in 
the variable £. This definition implies that a Gaussian state g is 
completely determined by the vector X of its first moments of 
the field operators, whose entries are given by Xj = Tr[gRj], 
and by the covariance matrix (CM) er, whose entries <Tjk are 
given by 

a jk = Tr [giRjRk + R k Rj)}/2 - X 3 X k . (3) 

Explicitly, the characteristic function of a Gaussian state 
with first moments X and CM <x is given by 

x (£) = c -èí T " T <TQí +í A' T r2£ _ (4) 

Gaussian states play a prominent role in practical realizations 
of continuous-variable (CV) quantum information protocols. 
They can be created and manipulated with relative ease with 
current technology fl·lll . and, thanks to their simple descrip- 
tion in terms of covariance matrices, provide a powerful and 
relevant theoretical framework for the investigation of funda- 
mental issues. 

All the unitary operations mapping Gaussian states into 
Gaussian states are generated by polynomials of the first and 
second order in the quadrature operators. First order opera- 
tions are just displacement operators Dç, which leave the CM 
unchanged while shifting the first moments. Such unitary op- 
erations, by which first moments can be arbitrarily adjusted, 
are manifestely local: this entails that first moments can play 
no role in the entanglement characterization of CV states and 
will be thus henceforth neglected, reducing the description of 
the states under exam to the CM er. On the other hand, uni- 
tary operations of the second order act, in Heisenberg picture, 
linearly on the vector R: R SR, where the matrix S sat- 
isfies S T ÍIS = 0. The set of such (real) matrices form the 
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real symplectic group Sp2n,n Íï^-Hóll . Therefore, these uni- 
tary operations are called symplectic operations. Symplectic 
operations act on a CM er by congruence: er i— > S 1 crS. 

Besides describing most unitary Gaussian operations cur- 
rently feasible in the experimental practice (namely beam- 
splitters, squeezers, and phase-shifters), the symplectic frame- 
work is fundamental in the theoretical analysis of CMs: for 
any physical CM er there exist a symplectic transformation 
S € Sp2n,n sucn that S t <tS = u, where 

v = ffi™=idiag(^-,^) . 

The quantities {vj }, uniquely determined for every CM er, are 
referred to as the symplectic eigenvalues of er, while v is said 
to be the Wüliamson normal form associated to er Od. It 
can be shown that, because of the canonical commutation rela- 
tions, the positivity of the density matrix g is equivalent to the 
following uncertainty relation for the symplectic eigenvalues 
of the CM describing a Gaussian state: 



Vj > 1 , for j = 1, . 



(5) 



The purity Tr [g 2 ] of a Gaussian state g with CM er and sym- 
plectic eigenvalues {vj} is simply givenby 



systems, as the action of partial transposition on CMs can be 
stated mathematically in very simple terms: the CM é> of the 
partially transposed state g with respect to, say, subsystem A, 
is simply obtained by switching the signs of the m momenta 
{pj} belonging to subsystem A ll22Íl : 



e> = TerT , with T = 



1 

-1 



(9) 



where la n stands for the 2ra-dimensional identity matrix. 
Even more remarkably, it has been proven that the PPT condi- 
tion is not only necessary, but as well sufficient for the sep- 
arability of (1 + n)-mode Gaussian states II2 2[ 12311 and of 
(m + ?i)-mode bisymmetric Gaussian states VM . thus pro- 
viding a powerful theoretical tool to detect quantum entangle- 
ment in these relevant classes of states. Let us notice that the 
(1 + n)-mode bipartitions encompass all the possible biparti- 
tions occurring in three-mode states. In analogy with Eq. (0, 
the PPT criterion can be explicitly expressed as a condition 
on the symplectic eigenvalues {pj] of the partially transposed 
CMe>: 



VA > 1 



for all 
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1. 



(10) 



Tr [g 2 ] = 1/VDct er = JJí 1 / 1 

3=1 



(6) 



The purity quantifies the degree of mixedness of the Gaussian 
state q, ranging from 1 for pure states to the limiting value 
for completely mixed states (due to the infinite dimension of 
the Hilbert space, no finite lower bound to the 2-norm of g ex- 
ist). Its conjugate S'l = 1 — Tr [g 2 ] is referred to as the linear 
entropy, ranging from for pure states to the limiting value 1 
for maximally mixed states. Another proper way of quantify- 
ing the mixedness of a state is provided by the von Neumann 
entropy Sv = — Tr [gln g]. The von Neumann entropy of a 
Gaussian state with CM er and symplectic eigenvalues {vj} 
reads lli 



3=1 



(7) 



with 



/(-) = ^lnf^-^lnf^l). (8) 



Let us now consider a (m + n)-mode bipartite Gaussian state 
i.e. a Gaussian state separated into a subsystem A of m modes, 
owned by party A, and a subsystem B of n modes, owned by 
party B. This state is associated to a 2(m + n)-dimensional 
CM er. Now, in general, for any bipartite quantum state g, the 
positivity of the partially transposed density matrix g, that is, 
the operator obtained from g by transposing the variables of 
only one of the two subsystems, is a necessary condition for 
the separability of the state. This condition thus goes under 
the name of "Positivity of Partial Transpose (PPT) criterion" 
1E0II2ÏI1 . This fact is especially useful when dealing with CV 



We finally mention that, in alternative to the PPT criterion, 
one can introduce an operational criterion based on a nonlin- 
ear map, that is independent of, and strictly stronger than the 
PPT condition ll25ll . In fact, this criterion is necessary and suf- 
ficient for separability of all (m + n)-mode Gaussian states of 
any m x n bipartions. 

For future convenience, let us define and write down the 
CM <J\... n of an n-mode Gaussian state in terms of two by 
two submatrices as 



/ -I 



Cl. ..71 = 



£12 



V 



^n—ln 
t n-ln °n ' 



(11) 



The symplectic eigenvalues v T of a two-mode CM er 12 are 
invariant under symplectic operations acting on er 12 . Start- 
ing from this observation, it has been shown that they can 
be retrieved from the knowledge of the symplectic invariants 
Det er 12 and A12 = Det <J\ + Det er 2 + 2Det £12, according 
to the following formula fÜ^Eóll : 



2/4 = A 12 T V A i2-4Detcr 12 . (12) 



The uncertainty relation Eq. l|5} imposes 
A12 -Dcter 12 < 1 . 



(13) 



Likewise, the symplectic eigenvalues v T of the CM ér 12 of 
the partially transposed state can be determined by partially 
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transposing such invariants and can thus be easily computed 

as 



2í>2 = Au =F 



Af 2 - 4Det cti2 



(14) 



where À12 = Dct <j\ + Dct er 2 — 2Dct £12. 
Let us finally observe that the quantities 



Al 



^Deto-j 



2 ^ Det i 

j<k 



are symplectic invariants for any number n of modes |2] 

We now move on to review in some detail the possible en- 
tanglement measures apt to quantify the entanglement of two- 
mode Gaussian states, upon which multipartite counterparts 
will be constructed in the following. 



B. Quantifying the entanglement of two-mode Gaussian states 

Thanks to the necessary and sufficient PPT criterion for 
separability, a proper measure of entanglement for two-mode 
Gaussian states is provided by the negativity Af, first intro- 
duced in Ref. 12811 . later thoroughly discussed and extended 
in Refs. I26l,l29ll to CV systems. The negativity of a quantum 
state q is defined as 



Af(g) 



(15) 



where g is the partially transposed density matrix and ||ò|| 1 = 
Tr|ó| stands for the trace norm of the hermitian operator d. 
This measures quantifies the extent to which g fails to be pos- 
itive. Strictly related to Af is the logarithmic negativity £jv, 
defined as £V = ln ||í?||i, which constitutes an upperbound 
to the distillable entanglement of the quantum state g and is 
related to the entanglement cost under PPT preserving opera- 
tions Both the negativity and the logarithmic negativity 
have been proven to be monotone under LOCC (local oper- 
ations and classical communication) l2ól I291 13111 . a crucial 
property for a bonafi.de measure of entanglement. Moreover, 
the logarithmic negativity possesses the agreeable property 
of being additive. For any two-mode Gaussian state g it is 
easy to show that both the negativity and the logarithmic neg- 
ativity are simple decreasing functions of the lowest symplec- 
tic eigenvalue í>_ of the CM of the partially transposed state 
MM: 



\g\\i = - = — => Aí(g) = max 



0, 



2£>_ 



Ej\f(o) = max [0, — lni/_] 



(16) 



(17) 



These expressions directly quantify the amount by which the 
necessary and sufficient PPT condition dlOl for separability is 
violated. The lowest symplectic eigenvalue V- of the partially 
transposed state cr thus completely qualifies and quantifies, in 



terms of negativities, the entanglement of a two-mode Gaus- 
sian state cr. For z>_ > 1 the state is separable, otherwise it 
is entangled; moreover, in the limit of vanishing z>_, the neg- 
ativities, and thus the entanglement, diverge. 

In the special instance of symmetric two-mode Gaussian 
states (i.e. of states with Det cr \ = Det <r 2 ), the entanglement 
of formation (EoF) Í32ll . can be computed as well l33ll . We 
recali that the EoF Ep of a quantum state g is defined as 

Ep(g) = min V Pl E(\^)) , (18) 

{pí,\4>í)} 

where E(\ipi}) denotes the von Neumann entropy Sv of the 
reduced density matrix of one party in the pure states E(\tpi)), 
namely the unique measure of bipartite entanglement for all 
pure quantum states (entropy of entanglement). The minimum 
in Eq. dl 8I > is taken over all the pure states realizations of g: 

Q = ^2Pi\^i)(ipi\ ■ 

i 

The asymptotic regularization of the entanglement of forma- 
tion coincides with the entanglement cost Ec(g), defined as 
the minimum number of singlets (maximally entangled anti- 
symmetric two-qubit states) which is needed to prepare the 
state g through LOCC Q. 

The optimal convex decomposition of Eq. dl 8i has been 
determined exactly for symmetric two-mode Gaussian states, 
and turns out to be Gaussian, that is, the absolute minimum 
is realized within the set of pure two-mode Gaussian states, 
yielding l33ll 



Ep = max [0, h{v-)} 



(19) 



with 



h(x) = 



{l + xf 



In 



■Ix 



(l-^) 2 
4x 



In 



(1-xf 



4.c 



(20) 

Such a quantity is, again, a monotonically decreasing function 
of £>_. Therefore it provides a quantification of the entangle- 
ment of symmetric states equivalent to the one provided by 
the negativities. This equivalence, regrettably, does not hold 
for general, mixed nonsymmetric states. In this case the EoF 
is not computable; nonetheless, it has been demonstrated that 
different entanglement measures induce different orderings of 
the states J35J]- This means that, depending on the measure of 
entanglement that one chooses, either the PPT-inspired nega- 
tivities or the entropy-based Gaussian measures (see below), a 
certain state can be more or less entangled than another given 
state. Clearly, this is neither a catastrophic nor an entirely un- 
expected result, but rather a consequence of the fact that, in 
general, for mixed states, different measures of entanglement 
may be associated to different conceptual and operational def- 
initions, and thus may measure different aspects of the quan- 
tum correlations present in a statistical mixture. 

In fact, restricting to the Gaussian framework, a special 
family of proper entanglement measures can be defined, shar- 
ing the agreeable property of being analytically computable 
in several instances of physical interest. The formalism of 
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Gaussian entanglement measures (Gaussian EMs), first intro- 
duced in Ref. 13611 . has been further developed and analysed 
tn Ref. 0. Such a formalism enables to define genèric Gaus- 
sian EMs of bipartite entanglement by applying the Gaussian 
convex roof, that is, the convex roof over pure Gaussian de- 
compositions only, to any bona fi.de measure of bipartite en- 
tanglement defined for pure Gaussian states. As already men- 
tioned, the optimization problem Eq. (II 81 for the computa- 
tion of the EoF of nonsymmetric two-mode Gaussian states 
has not yet been solved. However, the task can be somehow 
simplified by restricting to decompositions into pure Gaussian 
states only. The resulting measure, named "Gaussian EoF" in 
Ref. 13611 . is an upper bound to the true EoF and coincides 
with it for symmetric two-mode Gaussian states. 

In general, we can define a Gaussian EM Ge as follows. 
For any pure Gaussian state with CM er p , one has 

G e {<t p ) = E(\iP)) , (21) 

where E can be any proper measure of entanglement of pure 
states, defined as a monotonically increasing function of the 
entropy of entanglement (i.e. the von Neumann entropy of the 
reduced density matrix of one party). 

For any mixed Gaussian state g with CM er, one has ll36Íl 

G e (<t) = inf G e {<t p ) ■ (22) 

If the function E is taken to be exactly the entropy of entangle- 
ment, then the corresponding Gaussian EM defines the Gaus- 
sian entanglement of formation {Gaussian EoF) l3óll . Froman 
operational point of view, the Gaussian EoF is strictly related 
to the capacity of bosonic Gaussian channels Í37ÍI . Moreover, 
the Gaussian EoF is an entanglement monotone under Gaus- 
sian LOCC, a property that is shared by all Gaussian EMs 

S Hi- 
In general, the definition Eq. \22\ involves an optimization 
over all pure Gaussian states with CM er p smaller than the 
CM er of the mixed state whose entanglement one wishes to 
compute. This is a simpler optimization problem than that ap- 
pearing in the definition Eq. (II 81 of the true EoF, which, in 
CV systems, would imply considering decompositions over 
all, Gaussian and non-Gaussian pure states. Despite this sim- 
plification, in general the Gaussian EMs cannot be expressed 
in a simple closed form, even for two-mode Gaussian states. 
However, the Gaussian EMs have been computed analytically 
ll35ll for two relevant classes of, generally nonsymmetric, two- 
mode Gaussian states, namely the states of extremal - maxi- 
mal and minimal - negativity at fixed global and local puri- 
ties, referred to, respectively as Gaussian Maximally Entan- 
gled Mixed States (GMEMS) and Gaussian Least Entangled 
Mixed States (GLEMS) [kJ- In particular, the explicit ex- 
pression of the Gaussian EMs of the GLEMS will be crucial in 
the following because, as we are anout to show, any two-mode 
reduction of a three-mode pure Gaussian state is a GLEM. 

III. THREE-MODE GAUSSIAN STATES 

To begin with, let us set the notation and review the known 
results about three-mode Gaussian states of CV systems. We 



will refer to the three modes under exam as mode 1, 2 and 3. 
The two by two submatrices that form the CM er = er 12 3 of a 
three-mode Gaussian state are defined according to Eq. dl 1> . 
whereas the four by four CMs of the reduced two-mode Gaus- 
sian states of modes i and j will be denoted by (Tíj . Likewise, 
the local symplectic invariants Ay will be specified by the la- 
bels i and j of the modes they refer to, while, to avoid any 
confusion, the three-mode (global) symplectic invariant will 
be denoted by A = A123. Let us recali the uncertainty rela- 
tion Eq. dl 31 for two-mode Gaussian states: 

Ay - Det er.y < 1 . (23) 

As we have seen in the previous section, a complete qual- 
itative characterization of the entanglement of three-mode 
Gaussian state is possible because the PPT criterion is nec- 
essary and sufficient for their separability under any, partial 
or global, bipartition. This has lead to an exhaustive classi- 
fication of three-mode Gaussian states in five distinct classes 
fïlll . These classes take into account the fact that the modes 
1, 2 and 3 allow for three distinct biparti tions: 

• Class 1 : states not separable under all the three possi- 
ble bipartitions i x (jk) of the modes (fully inseparable 
states, possessing genuine multipartite entanglement). 

• Class 2: states separable under only one of the three 
possible bipartitions (one-mode biseparable states). 

• Class 3: states separable under only two of the three 
possible bipartitions (two-mode biseparable states). 

• Class 4: states separable under all the three possible 
bipartitions, but impossible to write as a convex sum of 
tripartite products of pure one-mode states (three-mode 
biseparable states). 

• Class 5: states that are separable under all the three pos- 
sible bipartitions, and can be written as a convex sum of 
tripartite products of pure one-mode states (fully sepa- 
rable states). 

Notice that classes 4 and 5 cannot be distinguished by par- 
tial transposition of any of the three modes (which is positive 
for both classes). States in class 4 stand therefore as nontriv- 
ial examples of tripartite entangled states of CV systems with 
positive partial transpose JTIll . It is well known that entangled 
states with positive partial transpose possess bound entangle- 
ment, that is, entanglement that cannot be distilled by means 
ofLOCC. 



A. Pure states 

We begin by focusing on pure three-mode Gaussian states, 
for which one has 

Dct er = 1 , A = 3 . (24) 

The purity constraint requires the local entropic measures of 
any 1 x 2-mode bipartitions to be equal: 

DeteTy = Det ert , (25) 
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with i, j and k different from each other. This general, well 
known property of the bipartitions of pure states may be easily 
proven resorting to the Schmidt decomposition. 

A first consequence of Eqs. (I24i and d25i is rather remark- 
able. Combining such equations one easily obtains 

(Ai 2 -Det<ri 2 ) + (A13 -Deto-is) 

+ (A 23 -Deto-as) =3 , (26) 

which, together with Inequality d23i . implies 

Ay = Dct + V i, ] : i£ ] . (27) 

The last equation shows that any reduced two-mode state of 
a pure three-mode Gaussian state saturates the partial uncer- 
tainty relation Eq. J23I . The states endowed with such a 
partial minimal uncertainty are states of minimal negativity 
for given global and local purities, Gaussian least entangled 
mixed states (GLEMS) 0IÏ3l- We reca U tnat by two-mode 
mixed Gaussian states of partial minimum Heisenberg uncer- 
tainty one means states that have one of the two symplectic 
eigenvalues equal to 1. States with both symplectic eigen- 
values equal to 1 are of course the pure Gaussian states of 
absolute minimum Heisenberg uncertainty. These definitions 
immediately extend to arbitary multimode Gaussian states. In 
this case, n-mode Gaussian states of m-partial minimum un- 
certainty are those that have m out of the n symplectic eigen- 
values equal to 1, with m < n. Notice that such a result could 
have also been inferred by invoking the reduction of (1 + n)- 
mode pure Gaussian states discussed in Ref. l38ll . first intro- 
duced in Ref. I39ÍI and proved at the covariance matrix level in 
Ref. lEoil . This implies that, through local unitàries (under any 
bipartition of the three modes), the state can be brought to the 
product of a two-mode squeezed state and of an uncorrelated 
vacuum. In turn, this implies that any of the three reduced 
two-mode CMs (resulting from the discarding of one mode) 
has one symplectic eigenvalue equal to 1 and is thus a GLEM. 

In fact, our simple proof, straightforwardly derived in terms 
of symplectic invariants, provides some further insight into 
the structure of CMs characterizing Gaussian states. What 
matters to our aims, is that the Standard form CM of Gaus- 
sian states is completely determined by their global and lo- 
cal invariants. Therefore, because of Eq. i25\ , the entangle- 
ment between any pair of modes embedded in a three-mode 
pure Gaussian state is fully determined by the local invariants 
Det ai , for / = 1,2,3, whatever proper measure we choose to 
quantify it l35ll . Furthermore, the entanglement of a (Ti\<Tjk 
bipartition of a pure three-mode state is determined by the en- 
tropy of one of the reduced states that is, once again, by the 
quantity Det er.j. Thus, the three local symplectic invariants 
Det <J\, Det er 2 and Det <x 3 fully determine the entanglement 
of any bipartition of a pure three-mode Gaussian state. We 
will show that they suffice to determine as well the genuine 
tripartite entanglement encoded in the state. 

For ease of notation, in the following we will denote by 
ai the local single-mode symplectic eigenvalues associated to 
mode l with CM af. 



Eq. i|6} shows that the quantities a; are simply related to the 
purities of the reduced single-mode states, the local purities 
pi, by the relation 

W = - • (29) 

ai 

Since the set {ai} fully determines the entanglement of any 
of the 1 x 2-mode and 1 x 1-mode bipartitions of the state, 
it is important to determine the range of the allowed vàlues 
for such quantities. This will provide a complete quantita- 
tive characterization of the entanglement of three-mode pure 
Gaussian states. To this aim, let us focus on the reduced two- 
mode CM <Ti2 and let us bring it (by local unitàries) in Stan- 
dard form I22L 14111 . so that Eq. (llit is recast in the form 

ai = diag{a ; , a ; }, l = 1,2; 
£12 = diag{ci2, di 2 } , (30) 

where C12 and du are the two-mode covariances, and, as 
we will show below, can be evaluated independently in pure 
three-mode Gaussian states. Notice that no generality is lost 
in assuming a Standard form CM, because the entanglement 
properties of any bipartition of the system are invariant under 
local (single-mode) symplectic operations. Now, Eqs. j25t 
and J24i may be recast as follows 

al = a\ + a\ + 2c 12 d 12 - 1 , (31) 
a 3 = (aia2 — c 12 )(aia2 — d 12 ) , (32) 

showing that we may eliminate one of the two covariances to 
find the expression of the remaining one only in terms of the 
three local inverse of the purities ai (mixednesses). Defining 
the quantity k as 

_ 1 + aj - a\- a\ 
K = c 12 d 12 = , (33) 

leads to the following condition on the covariance C\ 2 . 

4 \l i\2 i 2 2 2 21 2 1 2 n 

c 12 [{k - 1) + a x a 2 - % - a 2 \ c 12 + « = . 

a\a 2 

(34) 

Such a second order algebraic equation for cf 2 admits a posi- 
tive solution if and only if its discriminant S is positive: 

5 > . (35) 

After some àlgebra, one finds 

S = (ai + a 2 + a 3 + l)(ai + a 2 + a 3 - 1) 

x (ai + a 2 - a 3 + l)(ai - a 2 + a 3 + 1) 

x (— a\ + a 2 + a 3 + l)(ai + a 2 - a 3 - 1) 

x (ai — ü2 + a 3 — o-i + 0-2 +a>3 — 1) ■ (36) 

Aside from the existence of a real covariance C\ 2 , the fur- 
ther condition of positivity of <j\ 2 has to be fulfilled for a 
state to be physical. This amounts to impose the inequality 
a\a 2 — c\ 2 > 0, which can be explicitly written, after solving 
Eq. J34L as 



ai = v/Dcten 



(28) 



4 [2a\al - ((k - l) 2 + a\a\ -a\- aj)] > VS 
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This inequality is trivially satisfied when squared on both 
sides; therefore it reduces to 

2a\a\ - ((« - l) 2 + a\a\ - a\ - al) > . (37) 

Notice that conditions i35\ and d37i . although derived by 
assuming a specific bipartition of the three modes, are inde- 
pendent on the choice of the modes that enter in the consid- 
ered bipartition, because they are invariant under all possible 
permutations of the modes. Defining the parameters 

a\ = ai - 1 , (38) 

the Heisenberg uncertainty principle for single-mode states re- 
duces to 

a[>0 VI = 1,2,3. (39) 

This fact allows to greatly simplify the two previous condi- 
tions, which can be combined into the following triangular 
inequality 

\a'i -a' 3 \ < a' k < a\ + a' 3 . (40) 

Inequality (1401 is a condition invariant under all possible per- 
mutations of the mode indexes {i,j, k}, and, together with the 
positivity of each a[, fully characterizes the local symplectic 
eigenvalues of the CM of three-mode pure Gaussian states. It 
therefore provides a complete characterization of the entan- 
glement in such states. All Standard forms of pure three-mode 
Gaussian states and in particular, remarkably, all the possible 
vàlues of the logarithmic negativity between any pair of sub- 
systems, can be determined by letting a[, a' 2 and a' 3 vary in 
their range of allowed vàlues, as summarized in Fig. ^ 

Let us remark that Eq. J40i qualifies itself as an entropic 
inequality, as the quantities {a'j} are closely related to the pu- 
rities and to the von Neumann entropies of the single-mode 



reduced states. In particular the von Neumann entropies Syj 
of the reduced states are given by Syj = f(a'j + 1) = /(%), 
where the increasing convex entropic function f(x) has been 
defined in Eq. Now, Inequality J40i is strikingly anal- 
ogous to the well known triangle (Araki-Lieb) and subaddi- 
tivity inequalities for the von Neumann entropy (holding for 
general systems, see, e.g., fi]]), which in our case read 

l/K) " /K)l < /Ofc) < /(Oi) + /(%) • (41) 

However, as the different convexity properties of the functions 
involved suggest, Inequalities J40I and (14 1 i are not equiva- 
lent. Actually, as can be shown by exploiting the properties 
of the function f(x), the Inequalities d40i imply the Inequali- 
ties J4-1Í for both the leftmost and the rightmost parts. On the 
other hand, there exist vàlues of the local symplectic eigen- 
values {aj} for which Inequalities (14 1 i are satisfied but d40l > 
are violated. Therefore, the condition imposed by Eq. (I40> 
is stronger than the generally holding inequalities for the von 
Neumann entropy applied to pure states. 

Let us recali that the form of the CM of any Gaussian state 
can be simplified through local (unitary) symplectic opera- 
tions (that therefore do not affect the entanglement or mixed- 
ness properties of the state) belonging to Spf^,. Such reduc- 
tions of the CMs are called "Standard forms". For the sake of 
clarity, let us write the explicit Standard form CM of a genèric 
pure three-mode Gaussian state: 





«1 
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e + 

e 13 


o \ 







«1 





e l2 





e 13 




e 12 





«2 





e 23 















a 2 





e 23 









e 23 





«3 





V 





er 3 





e 23 





«3 / 



(42) 



with 



(at - aj) 2 - (a k - l) 2 



aj) - (a k + l) 2 



± 



(<z, + a,jf - (a k - l) 2 (a t + aj) - (a k + ïf 



4yc 



(43) 



By direct comparison with Eq. (67) in Ref. fïoll . it is imme- 
diate to verify that each two-mode reduced CM cr^ denotes 
a Standard form GLEMS with local purities /i; = a7/ and 
/ij = aj , and global purity ^ = [i k = . From our 
study it then turns out that, regarding the classification of Sec. 
mnuu, pure three-mode Gaussian states may belong either 
to class 5, in which case they reduce to the global three-mode 
vacuum, or to class 2, reducing to the uncorrelated product of 
a single-mode vacuum and of a two-mode squeezed state, or to 
class 1 (fully inseparable state). No two-mode or three-mode 
biseparable pure three-mode Gaussian states are allowed. 
Let us fïnally stress that, although useful in actual calcu- 



lations, the use of CMs in Standard form does not entail any 
loss of generality, because all the results derived in the present 
work do not depend on the choice of the specific form of the 
CMs, but only on invariant quantities, such as the global and 
local symplectic invariants. 



B. Mixed states 

The most general Standard form <r s f associated to the CM 
of any (generally mixed) three-mode Gaussian state can be 
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FIG. 1: (color online) Range of the entropic quantities a[ — /i; -1 — 1 for pure three-mode Gaussian states. The three parameters a' h with 
l = 1, 2, 3, have to vary inside the pyramid represented in plot (a) or, equivalently, for fixed vàlues of one of them, say aí, inside the shaded 
slice represented in plot (b), in order to determine the CM of a physical state, Eq. I42i . The expression of the boundary surfaces/curves come 
from the saturation of the triangular inequality <40t for all possible mode permutations. In particular, for the projected two-dimensional plot 
(b), the equations of the three boundaries are: I. a' 3 — a[ — a' 2 ; II. a' 3 — a[ + a' 2 \ III. a' 3 — a' 2 — a[. 



written as 



(o.x h h h\ 

oi h h 

h a 2 o h h 

h a 2 h h 

h h h ^3 

V/ 5 Ía fs h a 3 J 



(44) 



where the 12 parameters {cij} (inverse of the local purities) 
and {fj} (the covariances describing correlations between the 
modes) are only constrained by the Heisenberg uncertainty 
relations Eq. The possibility of this useful, general re- 
duction can be easily proven al ong the same lines as the two- 
mode Standard form reduction 14111 : by means of three local 
symplectic operations one can bring the three blocks a\, <r 2 
and <r 3 in Williamson form, thus making them insensitive to 
further local rotations (which are symplectic operations); ex- 
ploiting such rotations on mode 1 and 2 one can then diago- 
nalize the block £12 as allowed by its singular value decom- 
position; finally, one local rotation on mode 3 is left, by which 
one can cancel one entry of the block £13. Indeed, the re- 
sulting number of free parameters could have been inferred 
by subtracting the number of parameters of an element of 
Sp2,n © Sp2,n © Sp2,u (which is 9, as Sp2,B. has 3 inde- 
pendent generators) from the 21 entries of a genèric 6x6 
symmetric matrix. 



C. Symmetric states 

Among genèric Gaussian states, those endowed with some 
properties of symmetry under mode exchange play a special 
role for what concerns the structure of entanglement. In par- 
ticular, in a three-mode CV system, bisymmetric states are 



Gaussian states invariant under the exchange of two given 
modes (say 2 and 3) Q El ■ Their CM will be fhus of the 
form 




(45) 



Let mode 1 be entangled with the block of modes (23). It 
has been proven ll24ll42Íl fhat for such bisymmetric states the 
application of a local unitary (symplectic in phase space) op- 
eration on the block (23) concentrates the whole original mul- 
timode entanglement into the reduced state of a single pair of 
modes. Namely, in terms of the new modes {1,2', 3'}, the 
CM is transformed in a two-mode entangled state of modes 1 
and 2', tensor the uncorrelated single-mode state of mode 3', 
so that the original multimode entanglement can be quantified 
resorting to the well established theory of bipartite entangle- 
ment in two-mode Gaussian states USES Cl El 

The local symplectic transformation responsible for the 
unitary localization of the multimode entanglement is typi- 
cally realized by a simple beam splitter, if the CM is in Stan- 
dard form, with the single-mode blocks in their Williamson 
diagonal form. More generally, it may be a combination of 
beam splitters, phase shifters and squeezers. This type of 
entanglement localization is unitary and reversible, and thus 
completely different from the usual localization or concen- 
tration procedures that are based on measurements, as in the 
case of the "localizable entanglement" previously introduced 
for spin systems 1441 145ID . To reconstruct the original state, 
it suffices to let the discarded mode 3' interfere once more 
with mode 2' through the reversed beam splitter (that is, by 
applying the inverse symplectic operation). We remark that 
the unitary localizability is a property that extends to all 1 x n 
Gaussian states ll42ll . and to all m X n bisymmetric Gaussian 
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states |24], enabling two parties (owing two respective blocks 
of múltiple symmetric modes) to realize, by purely local con- 
trols, a perfect and reversible entanglement switch between 
two-mode and multimode quantum correlations. 

Three-mode Gaussian states which are invariant under the 
exchange of any two modes are said to be fully symmetric. 
They are trivially bisymmetric with respect to any 1x2 bipar- 
tition, meaning that each conceivable bipartite entanglement 
is locally equivalent to two-mode entanglement. In the Gaus- 
sian setting, these states are described by a CM ll24ll42Íl 



OL 


£ 


£ 


£ T 
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£ T 


£ T 
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where the local mixedness a = VDct a is the same for all the 
three modes. These states have been successfully produced in 
laboratory by quantum optical means iEóIEtIi . and exploited 
to implement quantum teleportation networks l48l Í49ll . Used 
as shared resources, they can be optimized with respect to lo- 
cal operations to realize CV teleportation with maximal non- 
classical fidel ity l50ll . quantum secret sharing IdÏI) . controlled 
dense coding l.52ll . and to solve CV Byzantine agreement Í53ll . 
Moreover, the structure of tripartite entanglement in this kind 
of states presents peculiar sharing properties íï^l . that are 
quite different from the properties of distributed entanglement 
among qubits and qudits [133, as will be discussed in detail in 
SecITcl 

We finally mention that the unitary localizability of entan- 
glement does not apply only to states with special symme- 
tries. For instance, for all pure three-mode Gaussian states, 
the 1x2 entanglement can be unitarily localized in any bi- 
partition. This fact holds f or g enèric pure Gaussian states of 
lxn biparti tions. QUlEj. 



entanglement (in particular, nonnegative on inseparable states 
and mono tònic under LOCC). 

It is natural to expect that Ineq. 1471 should hold for states 
of CV systems as well, despite the fact that they are defined 
on infinite-dimensional Hilbert spaces and can in principle 
achieve infinite entanglement, in particular the entanglement 
of distillation can become infinite in certain states of CV sys- 
tems; these states can be defined and constructed rigorously 
using the techniques of field theory and statistical mechan- 
ics for the description of systems of infinitely many degrees 
of freedom 15511 . In fact, one can show that the linearity of 
quantum mechanics, through the so-called no-cloning theo- 
rem ll56l l57l l58ll . prevents quantum correlations from being 
freely shareable, at striking variance with the behaviour of 
classical correlations fPïll . This entails that quantum entan- 
glement is "monogamous" ll59Íl . 

The crucial issue in contructing and proving the C V version 
of the CKW monogamy inequality is to find a proper measure 
of entanglement E, able to capture the trade-off between cou- 
plewise and tripartite correlations, quantitatively formalized 
by Ineq. (I47> . For qubit systems, such a measure is known 
as the tangle For Gaussian states of CV systems, this 
problem has been recently solved in Ref. fÏ2ll . where the CV 
analogue of the tangle has been defined and exploited to ob- 
tain a proof of the monogamy inequality J47i for all Gaussian 
states of three modes, and for all symmetric Gaussian states 
of systems with an arbitrary number of modes. Following the 
approach of Ref. ÍT3l . we recali now the notation leading to 
the definition of the continous-variable tangle, and provide a 
detailed proof of the CKW monogamy inequality obeyed by 
all three-mode Gaussian states. 



B. The continuous-variable tangle 



IV. GENUINE TRIPARTITE ENTANGLEMENT AND 
ENTANGLEMENT SHARING 

In this section we approach in a systematic way the question 
of distributing quantum correlations among three parties glob- 
ally prepared in a (pure or mixed) three-mode Gaussian state, 
and we deal with the related problem of quantifying genuine 
tripartite entanglement in such a state. 

A. Entanglement sharing 

The key ingredient of our analysis is the so-called sharing 
or monogamy inequality, first introduced by Coffman, Kundu, 
and Wootters (CKW) |J] for systems of three qubits, and re- 
cently extended to systems of n qubits by Osborne and Ver- 
straete J4J. The CKW monogamy inequality for a three-party 
system can be written as follows: 

£i|0'*) _ E l \ j - E^ k > , (47) 

where i, j, k denote the three elementary parties (modes in 
a CV system), and E refers to a proper measure of bipartite 



The continuous-variable tangle E T is formally defined as 
follows 1121 . For a genèric pure state of a (1 + 7V)-mode 
CV system, one has 

£ T W0=ln 2 ||p||i, P=\Í>M- (48) 

This is a proper measure of bipartite entanglement, being a 
convex, increasing function of the logarithmic negativity E/sí, 
equivalent to the entropy of entanglement on pure states. For 
a pure Gaussian state |^>) with CM a p , it is easy to find that 



E T (a p ) = arcsinh 2 ^ ^ , (49) 



where /ii = 1/ \/Dct <J\ is the local purity of the reduced state 
of mode 1, described by a CM ct\ (we are considering a most 
general lxn bipartition). Def. ( 148 \ is naturally extended to 
genèric mixed states p of (n + l)-mode CV systems through 
the convex-roof formalism lóÒTl . Namely, 

E T (p) = inf Y^piEr^i), (50) 
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where the infimum is taken over all convex decompositions 
of p in terms of pure states {IV'i)}- If the index i is contin- 
uous, the sum in Eq. J50l > is replaced by an integral, and the 
probabilities {pi} by a probability distribution ir(ip). 

Next, it is important to recali that for two qubits the tangle 
can be equivalently defined as the convex roof of the squared 
negativity lóHl . because the latter coincides with the concur- 
rence for pure two-qubit states Ió3l . Then, Eq. J50l l states 
that the convex roof of the squared logarithmic negativity de- 
fines the proper continuous-variable tangle, or, in short, the 
contangle E T (p) fï^l . One could have defined the contan- 
gle using the convex roof extension of the squared negativ- 
ity as well. The two defini tions are, in fact, equivalent to 
the aim of quantifying distributed entanglement, because the 
squared negativity is a convex function of the squared loga- 
rithmic negativity I3l 1 1 311 - The nice feature of using specif- 
ically the squared logarithmic negativity lies in the fact that 
from a computational point of view the logarithm accounts in 
a straightforard way for the infinite dimensionality of the un- 
derlying Hilbert sapce fÏ2Íl . We will prové in the following 
that the contangle satisfies the CKW monogamy inequality 
for all three-mode Gaussian states. Viceversa, one can easily 
show that any continuous-variable tangle defined in terms of 
the (not squared) negativity or of the entanglement of forma- 
tion fails to satisfy the CKW monogamy inequality in general 
ffHl . This situation is to some extent reminiscent of the case 
of qubit systems, for which the CKW monogamy inequal- 
ity holds using the tangle, defined as the convex roof of the 
squared concurrence or of the squared negativity lólll . but 
fails if one chooses alternative defini tions based on the convex 
roof of other equivalent measures of bipartite entanglement, 
such as the concurrence itself or the entanglement of forma- 
tion H. 

From now on, we restrict our attention to Gaussian states. 
Any multimode mixed Gaussian state with CM cr, admits a 
decomposition in terms of pure Gaussian states only. The in- 
fimum of the average contangle, taken over all pure Gaussian 
state decompositions, defines the Gaussian contangle G T 

GJo) = inf / n(dcr p )E T (a p ) . (51) 

{7r(ti<TP),<TP} J 

It follows from the convex roof construction that the Gaus- 
sian contangle G T {cr) is an upperbound to the true contangle 
E T (cr) (because the latter can be in principle minimized over 
a non-Gaussian decomposition): 

E T {a) < G r (<x) , (52) 

and it can be shown that G T {cr) is a bipartite entanglement 
monotone under Gaussian local operations and classical Com- 
munications (GLOCC) I35II3IÜI . The Gaussian contangle can 
be expressed in terms of CMs as 

G r (cr) = inf E T {(T P ) , (53) 

where the infimum runs over all pure Gaussian states with CM 
cr p < cr. Let us remark that, if cr s denotes a mixed symmetric 
(lx l)-mode Gaussian state, then the decomposition of cr s in 



terms of an ensemble of pure Gaussian states is the optimal 
one l33ll . which means that the Gaussian contangle coincides 
with the true contangle. Moreover, the optimal pure-state CM 
cr p minimizing G T (<r s ) in Eq. d53l > is characterized by having 
v-{cr p s ) = ü-(a s ) I33H36I1 . The fact that the smallest sym- 
plectic eigenvalue is the same for both partially transposed 
CMs entails for symmetric two-mode Gaussian states that 

E T {cr s ) = G T {a s ) = [max{0,-mi>_(<T s )}] 2 ■ (54) 

Finally, of course E T = G T as well in all pure Gaussian states 
of 1 x n bipartitions. 

C. Monogamy inequality for all three-mode Gaussian states 

We now provide the detailed proof, first derived, among 
other results, in Ref. 112Í1 . that all three-mode Gaussian states 
satisfy the CKW monogamy inequality ( I47i . using the (Gaus- 
sian) contangle to quantify bipartite entanglement. The inter- 
mediate steps of the proof will be then useful for the subse- 
quent computation of the residual genuine tripartite entangle- 
ment, as we will show in Sec. I1V Dl 

We start by considering pure three-mode Gaussian states, 
whose Standard form CM cr p is given by Eq. J42l >. As dis- 
cussed in Sec. Ull Al all the properties of entanglement in 
pure three-mode Gaussian states are completely determined 
by the three local purities. Reminding that the mixednesses 
ai = have to vary constrained by the triangle inequality 
040b . in order for a v to represent a physical state, one has 

\a 3 - Ofe| + 1 < ai < dj + ük — 1 ■ (55) 

For ease of notation let us rename the mode indices so that 
{i,j, k} = {1, 2, 3}. Without any loss of generality, we can 
assume a\ > 1. In fact, if a\ = 1 the first mode is not cor- 
related with the other two and all the terms in Ineq. \A1\ are 
trivially zero. Moreover, we can restrict the discussion to the 
case of both the reduced two-mode states er 12 and <x 13 being 
entangled. In fact, if e.g. <r 13 denotes a separable state, then 

i?r' 2 < Er^ 23 ^ because tracing out mode 3 is a LOCC, and 
thus the sharing inequality is automatically satisfied. We will 
now prové Ineq. d47l in general by using the Gaussian contan- 
gle, as this will immediately imply the inequality for the true 
contangle as well. In fact, G} I(23) (itP) = £^ l(23) (a p ), but 

Gl ll (cr) > E$\tr), 1 = 2,3. 

Let us proceed by keeping a\ fixed. From Eq. \A9\ . it fol- 
lows that the entanglement between mode 1 and the remaining 
modes, £^ 23 ' = arcsinb 2 ^ a\ — 1, is constant. We must 
now prové that the maximum value of the sum of the 1 1 2 and 
1 1 3 bipartite entanglements can never exceed E\ "- 23 - 1 , at fixed 
local mixedness a\. Namely, 

max Q < arcsinb 2 \/a 2 — 1 , (56) 

s,d 

where a = ai (from now on we drop the subscript "1"), and 
we have defined 

Q = G]} 2 (<t p ) + G]} 3 (<t p ) . (57) 
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The maximum in Eq. (I56> is taken with respect to the "center 
of mass" and "relative" variables s and d that replace the local 
mixednesses a-i and 03 according to 



d = 



a-2 + Q3 
2 

a.2 - a-3 



(58) 



(59) 



The two parameters s and d are constrained to vary in the 
region 



s 5 



a+l 



\d\ < 



a 2 -l 
As 



(60) 



Ineq. (I60i combines the triangle inequality J55i with the con- 
dition of inseparability for the states of the reduced biparti- 
tions (1|2) and (1|3) 0. 



We recali now, as stated in Sec. MI Al that each cru, 
l = 2, 3, is a state of partial minimum uncertainty (GLEMS 
fïcln . For this class of states the Gaussian measures of en- 
tanglement, including G T , can be computed explicitely 
yielding 



Q = arcsinh 



\J m 2 (a, s, d) — 1 



arcsinh^ 



\J m 2 (a, s, —d) — 1 



(61) 



where m = m_ if D < 0, and m = m+ otherwise (one has m+ = to_ for D = 0). Here: 



(s -d) 2 -l 



2a 2 (l + 2.s 2 + 2d?) - (4s 2 - l)(4d 2 - 1) - a 4 - y/ó 
4(s - d) 



L> = 2(s -d)-J2 [k 2 _ + 2k + + |fc_|(fc 2 + 8fc+) 1 /2] , 

fc± = a 2 ±{s + d) 2 , (62) 

and the quantity 

5={a-2d - l)(a -2d+ \){a + 2d- l)(a + 2d + l)(a - 2s - l)(a - 2s + l)(a + 2s - l)(a + 2s + 1) 



is the same as in Eq. J35I . Note (we omitted the explicit de- 
pendence for brevity) that each quantity in Eq. i62\ is a func- 
tion of (a,s,d). Therefore, to evaluate the second term in 
Eq. ( 16 U each d in Eq. i62\ must be replaced by —d. 

Studying the derivative of m T with respect to s, it is analyt- 
ically proven that, in the whole range of parameters {a, s, d} 
defined by Ineq. J60l >. both m_ and m + are monotonically 
decreasing functions of s. The quantity Q is then maximized 
over s for the limiting value 



(63) 



r 



<Ti3, the Gaussian contangle is defined in terms of m_. The 
latter, in turn, acquires the simple form 



m_(a,s min ,d) = 



1 + 3a + 2d 
3 + a - 2d 



(64) 



Consequently, the quantity Q turns out to be an even and con- 
vex function of d, and this fact entails that it is globally maxi- 
mized at the boundary 



| d | = rf max = 



a-1 



(65) 



This value of s corresponds to three-mode pure Gaussian 
states in which the state of the reduced bipartition 2|3 is al- 
ways separable, as one should expect because the bipartite 
entanglement is maximally concentrated in the states of the 
1|2 and 1 1 3 reduced bipartitions. With the position Eq. ( I63l >. 
the quantity D defined in Eq. Í621 can be easily shown to be 
always negative. Therefore, for both reduced CMs er 12 and 



We finally have that 

gmax 



Q[a,s = s min ,d-- 
arcsinh 2 \J a 2 — 1 



±d n 



(66) 



which implies that in this case the sharing inequality J47t is 
exactly saturated and the genuine tripartite entanglement is 
consequently zero. In fact this case yields states with «2 = ai 
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and 0,3 = 1 (if d = <i max ), or 03 = a\ and 02 = 1 (if 
d = — <i max ), z.e. tensor products of a two-mode squeezed 
state and a single-mode uncorrelated vacuum. Being Ç; max 
from Eq. i66\ the global maximum of Q, Ineq. (I56> holds true 
and the monogamy inequality (I47> is thus proven for any pure 
three-mode Gaussian state, choosing either the Gaussian con- 
tangle G T or the true contangle E T as measures of bipartite 
entanglement JÏ2Í1 . 

The proof immediately extends to all mixed three-mode 
Gaussian states <x, but only if the bipartite entanglement is 
measured by G T (cr) Í63ll . Let {Tr(da-f n ) 7 <J v m \ be the ensem- 
ble of pure Gaussian states minimizing the Gaussian convex 
roof in Eq. (15 li : then, we have 



G^(a) = J .(K)Gj (,t) K) 



> J .KJIG^'KJ + Gf «)] (67) 

where we exploited the fact that the Gaussian contangle is 
convex by construction. This concludes the proof of the CKW 
monogamy inequality J47i for all three-mode Gaussian states. 

We close this subsection by discussing whether the CKW 
monogamy inequality can be generalized to all Gaussian 
states of systems with an arbitrary number n + 1 of modes. 
Namely, we want to prové that 



E 1 ^ > 



(68) 



z=i 



Establishing this result in general is a highly nontrivial task, 
but it can be readily proven for all symmetric multimode Gaus- 
sian states JÏ2Í1 . In a fully symmetric n + 1-mode Gaussian 
state all the local purities are degenerate and reduce to a single 
parameter ai oc : 



.12 



"lo 



(69) 



As in the three-mode case, due to the convexity of G T , it will 
suffice to prové Eq. (16 8 1 for pure states, for which the Gaus- 
sian contangle coincides with the true contangle in every bi- 
partition. For any n and for a\ oc > 1 (for ai oc = 1 we have a 
product state), one has that 



E Wu-,j n ) =l n 2 (afoc _ ^Jal c -1) (70) 
is independent of n, while the total two-mode contangle 



nEÍ\ jl 
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o>L{n + l)-l 



l)(a 2 oc (n + l)2-(n-l)2) /„ 



"loc 



(71) 



is a monotonically decreasing function of the integer n at fixed 
a; oc . Because the sharing inequality trivially holds for n = 1, 
it is inductively proven for any n. This result, together with 



extensive numerical evidence obtained for randomly gener- 
ated non-symmetric 4-mode Gaussian states fÏ2ll . strongly 
supports the conjecture that the CKW monogamy inequality 
holds true for all multimode Gaussian state, using the (Gaus- 
sian) contangle as a measure of bipartite entanglement. How- 
ever, at present, a fully analytical proof of this conjecture is 
still lacking. 



D. Residual contangle, genuine tripartite entanglement, and 
monotonicity 

The sharing constraint leads naturally to the definition of 
the residual contangle as a quantifier of genuine tripartite en- 
tanglement (arravogliament) in three-mode Gaussian states, 
much in the same way as in systems of three qubits yj . How- 
ever, at variance with the three-qubit case, here the residual 
contangle is partition-dependent according to the choice of 
the reference mode, with the exception of the fully symmetric 
states. A bonafi.de quantification of tripartite entanglement is 
then provided by the minimum residual contangle 11211 



E^ k = min ES (jk) - ES j - Ef 



(72) 



where the symbol (i, j, k) denotes all the permutations of the 

three mode indexes. This definition ensures that Er^ k is in- 
variant under all permutations of the modes and is thus a gen- 
uine three-way property of any three-mode Gaussian state. We 
can adopt an analogous definition for the minimum residual 
Gaussian contangle G r T es (see Fig. |2]for a pictorial represen- 
tation): 



G res = Qi\j\k = min 
(i,j,k) 



(73) 



One can verify that 



i\m 



Gf) - (G J T 



■j\m 



G ] T lk ) > (74) 



if and only if > üj, and therefore the absolute minimum 
in Eq. d72t is attained by the decomposition realized with re- 
spect to the reference mode l of smallest local mixedness ai, 
i.e. for the single-mode reduced state with CM of smallest 
determinant. 

The residual (Gaussian) contangle must be nonincreasing 
under (Gaussian) LOCC in order to be a proper measure of 
tripartite entanglement. The monotonicity of the residual tan- 
gle was proven for three-qubit pure states in Ref. l64ll . In the 
CV setting, it has been shown in Ref. lfÏ2l that for pure three- 
mode Gaussian states the residual Gaussian contangle Eq. d73l 
is an entanglement monotone under tripartite GLOCC, and 
that it is nonincreasing even under probabilistic operations, 
which is a stronger property than being only monotone on 
average. Therefore the Gaussian contangle G r T es defines (to 
the best of our knowledge) thefirst measure, proper and com- 
putable, of genuine multipartite (specifically, tripartite) entan- 
glement in Gaussian states of CV systems. It is worth noting 
that the minimum in Eq. d73i . that at first sight might appear 
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FIG. 2: (color online) Pictorial representation of Eq. 1731 . defin- 
ing the residual Gaussian contangle G% £a of genèric (nonsymmetric) 
three-mode Gaussian states. G r T es quantifies the genuine tripartite 
entanglement shared among mode 1 (•), mode 2 (■), and mode 3 (,*.). 
The optimal decomposition that realizes the minimum in Eq. {73) is 
always the one for which the CM of the reduced state of the reference 
mode has the smallest determinant. 



a redundant requirement, is physically meaningful and math- 
ematically necessary. In fact, if one chooses to fix a reference 
partition, or to take e.g. the maximum (and not the minimum) 
over all possible mode permutations in Eq. (I73> . the resulting 
"measure" is not monotone under GLOCC and thus is defi- 
nitely not a measure of tripartite entanglement. 

We now work out in detail an explicit application, by de- 
scribing the complete procedure to determine the genuine tri- 
partite entanglement in a pure three-mode Gaussian state cr p . 

(i) Determine the local purities. The state is globally pure 

(Dct <j p = 1); therefore, the only quantities needed 
for the computation of the tripartite entanglement are 
the three local mixednesses au defined by Eq. d28i . 
of the single-mode reduced states <Ti, l = 1, 2, 3 (see 
Eq. (O)- Notice that the global CM a p needs not to 
be in the Standard form (142 > . as the single-mode de- 
terminants are local symplectic invariants fiDIl , From 
an experimental point of view, the parameters ai can 
be extracted from the CM using the homodyne tomo- 
graphic reconstruction of the state l65ll : or they can be 
directly measured with the aid of single photon detec- 
tors (óHEzTi 

(ii) Find the minimum. From Eq. (I74> . the minimum in the 

definition (I73> of the residual Gaussian contangle G r T es 
is attained in the partition where the bipartite entangle- 
ments are decomposed choosing as reference mode l the 
one in the single-mode reduced state of smallest local 
mixedness ai = a m i n . 

(iii) Check range and compute. Given the mode with small- 

est local mixedness a m i n (say, for instance, mode 1) 
and the parameters s and d defined in Eqs. A58I59I . if 
a-min = 1 then mode 1 is uncorrelated from the others: 
Qres _ q jf^ instem a m i n > 1 then 




FIG 3: (color online) Three-dimensional plot of the residual Gaus- 
sian contangle G r T es (cr p ) in pure three-mode Gaussian states cr p , de- 
termined by the three local mixedness ai, l — 1,2, 3. One of the 
local mixedness is kept fixed (a\ — 2). The remaining ones vary 
constrained by the triangle inequality 1551 . as depicted in Fig. Qb). 
The explicit expression of G r T es is given by Eq. I27i . See text for 
further details. 
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l)/4s then Gr 3 



0. Otherwise, all terms 



G r T es {(T p ) = arcsinh' 



2 

min 



in G r T es Eq. M2>\ are nonvanishing. 

The residual Gaussian contangle Eq. d73i in genèric pure 
three-mode Gaussian states is plotted in Fig. [3]as a function of 
0,2 and a 3 , at constant ai = 2. For fixed a\, it is interesting to 
notice that G r T es is maximal for 02 = 03, i.e. for bisymmetric 
states. Notice also how the residual Gaussian contangle of 
these bisymmetric pure states has a cusp for a\ = 0.2 = 03. 
In fact, from Eq. d74> . for 02 = 03 < ai the minimum in 
Eq. d73l is attained decomposing with respect to one of the 
two modes 2 or 3 (the result is the same by symmetry), while 
for ct2 = 03 > ai mode 1 becomes the reference mode. 

For genèric mixed three-mode Gaussian states, a quite cum- 
bersome analytical expression for the 1|2 and 1 1 3 Gaussian 
contangles may be written l35t l36ll . involving the roots of a 
fourth order polynomial, but the optimization appearing in the 
computation of the 1 1 (23) bipartite Gaussian contangle (see 
Eq. ( 153» has to be solved only numerically. However, ex- 
ploiting techniques like the unitary localization l24ll described 
in Sec. Ull Cl and results like that of Eq. \5A\ , closed expres- 
sions for the residual Gaussian contangle can be found as well 
in relevant classes of mixed three-mode Gaussian states en- 
dowed with some symmetry constraints. Interesting examples 
of these states and the investigation of their physical proper- 
ties will be discussed in Sec. [V] 

As an additional remark, let us recali that, although the en- 
tanglement of Gaussian states is always distillable with re- 
spect to 1 x N biparti tions 12311 . they can exhibit bound entan- 
glement in 1 x 1 x 1 tripartitions lïHl . In this case, the resid- 
l)/4s then Gt ' 2 = 0. Instead, if ual contangle cannot detect tripartite PPT entangled states. 



1 



d), (75) 



with Q = Gr' 2 + Gj.' 3 defined by Eqs. doTl62l . Note 
that if d < — (a 



14 



For example, the residual contangle in three-mode bisepara- 
ble Gaussian states (class 4 of Ref. fiUp is always zero, be- 
cause those bound entangled states are separable with respect 
to all (1 x 2)-modebipartitions. In this sense we can correctly 
regard the residual contangle as an estimator of distillable tri- 
partite entanglement in fully inseparable three-mode Gaussian 
states. However, we remind that this entanglement can be dis- 
tilled only resorting to non-Gaussian LOCC 1 68], since dis- 
tilling Gaussian states with Gaussian operations is impossible 

MfmïïÈ- 



V. SHARING STRUCTURE OF TRIPARTITE 
ENTANGLEMENT 

We are now in the position to analyze the sharing structure 
of CV entanglement in three-mode Gaussian states by taking 
the residual Gaussian contangle as a measure of tripartite en- 
tanglement, in analogy with the study done for three qubits 
116411 using the residual tangle yfl. 

The first task we face is that of identifying the three-mode 
anàlogues of the two inequivalent classes of fully inseparable 
three-qubit states, the GHZ state fïH 



IV'GHZ) = -J=(|000> + |111» 



and the W state H 
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(|001) + |010) + |100>) 



(76) 



(77) 



These states are both pure and fully symmetric, i.e. invariant 
under the exchange of any two qubits. On the one hand, the 
GHZ state possesses maximal tripartite entanglement, quanti- 
fied by the residual tangle I3l· l64ll . with zero couplewise entan- 
glement in any reduced state of two qubits reductions. There- 
fore its entanglement is very fragile against the loss of one or 
more subsystems. On the other hand, the W state contains the 
maximal two-party entanglement in any reduced state of two 
qubits [64] and is thus maximally robust against decoherence, 
while its tripartite residual tangle vanishes. 



A. CV GHZ/W states 

To define the CV counterparts of the three-qubit states 
IV'Ghz) and \ipw}> one must start from the fully symmetric 
three-mode CM cr s of Eq. J46i . Surprisingly enough, in sym- 
metric three-mode Gaussian states, if one aims at maximizing, 
at given single-mode mixedness a; oc = a, either the bipartite 
entanglement G % P in any two-mode reduced state (i.e. aiming 
at the CV WMike state), or the genuine tripartite entanglement 
G r T es {i.e. aiming at the CV GHZ-like state), one finds the 
same, unique family of pure symmetric three-mode squeezed 
states q-p. These states, previously known as CV "GHZ-type" 
states Í48T. 1731 F74I1 . can be indeed defined for genèric n-mode 
systems. They constitute an ideal test-ground for the study of 
the scaling of multimode CV entanglement with the number 



of modes. This analysis can be carried out via nested appli- 
cations of the procedure of unitary localization Í24L l4*3l . re- 
viewed in Sec. IIII Cl For systems of three modes, they are 
described by a CM crf of the form Eq. J46i . with a = al2, 
e = diag{e + , e~} and I42Í1 



e± = 



1 ± yV - 1) (9a 2 - 1) 
4a 



(78) 



ensuring the global purity of the state. For self-explaining 
reasons, we choose to name these states "CV GHZ/W states" 
JÏ2ÍI . and denote their CM by crf HZ ^ W . In the limit of infi- 
nite squeezing (a — ► 00), the CV GHZ/W state approaches 
the proper (unnormalizable) continuous-variable GHZ state 
/ dx\x 1 x, x), a simultaneous eigenstate of total momentum 
P1+P2+P3 andof all relative positions — ij = 1,2,3), 
with zero eigenvalues l75ll . 

The residual Gaussian contangle of GHZ/W states of finite 
squeezing takes the simple form lfÏ2ll 



G r T es {(rf HZ/w ) = arcsinh 2 y/a 2 - 1 



l· l 2 



3a 2 - 1 - V9a 4 - 10a 2 + 1 



(79) 



It is straightforward to see that G r T es {<T v s ) is nonvanishing as 
soon as a > 1. Therefore, the GHZ/W states belon g to 
the class of fully inseparable three-mode states fill l42t I48L 
0IÏÍ. We finally recali that in a GHZ/W state the resid- 
ual Gaussian contangle G r T es Eq. M3\ coincides with the true 
residual contangle E r T es Eq. J72i . This property clearly holds 
because the Gaussian pure-state decomposition is the optimal 
one in every bipartition, due to the fact that the global three- 
mode state is pure and the reduced two-mode states are sym- 
metric. 



B. T states 

The peculiar nature of entanglement sharing in C V GHZ/W 
states is further confirmed by the following observation. If 
one requires maximization of the 1x2 bipartite Gaussian 
contangle Gr under the constraint of separability of all 
the reduced two-mode states, one finds a class of symmet- 
ric mixed states characterized by being three-mode Gaussian 
states of partial minimum uncertainty. They are in fact char- 
acterized by having their smallest symplectic eigenvalue equal 
to 1, and represent thus the three-mode generalization of two- 
mode symmetric GLEMS H fïoll3l . 

We will name these states T states, with T standing for 
tripartite entanglement only. They are described by a CM crj 
of the form Eq. J46l >. with a = ala, £ = diag{e + , e~} and 



5 + y/9a 2 (a 2 - 2) + 25 
4a 



5 - 9a 2 + ^a 2 (a 2 - 2) + 25 



12a 



(80) 
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The T states, like the GHZ/W states, are determined only 
by the local mixedness a, are fully separable for a = 1, and 
fully inseparable for a > 1. The residual Gaussian contangle 
Eq. (I73> can be analytically computed for these mixed states 
as a function of a. First of all one notices that, due to the com- 
plete symmetry of the state, each mode can be chosen indiffer- 
ently to be the reference one in Eq. ( I73> . Being the lxl entan- 
glements all zeroby construction, G r T es = Gr ■ The 1x2 
bipartite Gaussian contangle can be in turn obtained exploit- 
ing the unitary localization procedure (see Sec. IIII C> . Let us 
choose mode 1 as the reference mode and combine modes 2 
and 3 at a 50:50 beam splitter, a local unitary operation with 
respect to the bipartition 1 1 (23) that defines the transformed 
modes 2' and 3'. The CM <j t s of the state of modes 1, 2', and 
3' is then written in the following block form: 




(81) 



where mode 3' is now disentangled from the others. Thus 



G^ 23 \c 



(82) 



Moreover, the reduced CM cr^ of modes 1 and 2' defines a 
nonsymmetric GLEM BO with 



Dct <J\ 

Dct er 2 



2)a 2 



Dct cria' = — ( 3a 



25-1 
3 



i (3a 2 + ^W 2 

(3a 2 - v/9 (a 2 - 2) a 2 + 25 + 3) 



and it has been shown that the Gaussian contangle is com- 
putable in two-mode GLEMS l35ll . After some àlgebra, one 
finds the complete expression of G r T es for T states: 





FIG. 4: (color online) Plot, as a function of the single-mode mixed- 
ness a, of the tripartite residual Gaussian contangle G r T es Eq. 1791 
in the CV GHZ/W states (dashed line); in the T states Eq. 483} 
(sòlid line); and in 50000 randomly generated mixed symmetric 
three-mode Gaussian states of the form Eq. 1461 (dots). The GHZ/W 
states, that maximize any bipartite entanglement, also achieve max- 
imal genuine tripartite quantum correlations, showing that CV en- 
tanglement distributes in a promiscuous way in symmetric Gaussian 
states. Notice also how all random mixed states have a nonnegative 
residual Gaussian contangle. This confirms the results presented in 
Ref. and discussed in detail and extended in Sec. lIVCl on the 
strict validity of the CKW monogamy inequality for CV entangle- 
ment in three-mode Gaussian states. 



restricting to pure Gaussian decompositions only in the defi- 
nition Eq. i73\ of the residual Gaussian contangle. In fact, for 
T states the relation G r T es (cr^) > £Ç es (of ) holds due to the 
symmetry of the reduced two-mode states, and to the fact that 
the unitarily transformed state of modes 1 and 2' is mixed and 
nonsymmetric. The crucial consequences of this result for the 
structure of the entanglement trade-off in Gaussian states will 
be discussed further in the next subsection. 



G*; es (cr^) = arcsinh 



25i? - 9a 4 + 3Ra 2 + 6a 2 - 109 
(81a 8 - 432a 6 + 954a 4 - 1704a 2 + 2125 
(3a 2 - 11) (3a 2 - 7) (3a 2 + 5) i?) 2 V2 
[18 (3a 2 - R + 3)}~^ i , (83) 



with R = v/9a 2 (a 2 - 2) + 25. 

What is remarkable about T states is that their tripartite 
Gaussian contangle Eq. J83I is strictly smaller than the one 
of the GHZ/W states Eq. 03 for any fixed value of the local 
mixedness a, that is, for any fixed value of the only parameter 
(operationally related to the squeezing of each single mode) 
that completely determines the CMs of both families of states 
up to local unitary operations. This hyerardical behavior of 
the residual contangle in the two classes of states is illustrated 
in Fig. [4] Notice that this result cannot be an artifact caused by 



C. Promiscuous continuous-variable entanglement sharing 

The above results, pictorially illustrated in Fig. [4] lead to 
the conclusion that in symmetric three-mode Gaussian states, 
when there is no bipartite entanglement in the two-mode re- 
duced states (like in T states) the genuine tripartite entangle- 
ment is not enhanced, but frustrated. More than that, if there 
are maximal quantum correlations in a three-party relation, 
like in GHZ/W states, then the two-mode reduced states of 
any pair of modes are maximally entangled mixed states. 

These findings, unveiling a major difference between 
discrete-variable (mainly qubits) and continuous-variable sys- 
tems, establish the promiscuous nature of CV entanglement 
sharing in symmetric Gaussian states 1IT2I [I3I1 . Being as- 
sociated with degrees of freedom with continuous spectra, 
states of CV systems need not saturate the CKW inequality 
to achieve maximum couplewise correlations. In fact, without 
violating the monogamy constraint Ineq. ( I47K pure symmetric 
three-mode Gaussian states are maximally three-way entan- 
gled and, at the same time, maximally robust against the loss 
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of one of the modes. This preselects GHZ/W states also as 
optimal candidates for carrying quantum information through 
a lossy channel, being, for their intrinsic entanglement struc- 
ture, less sensitive to decoherence effects, as we will show in 
Sec. EU 

As an additional remark, let us mention that, qui te naturally, 
not all three-mode Gaussian states (in particular nonsymmet- 
ric states) are expected to exhibit a promiscuous entanglement 
sharing. Further investigations to clarify the sharing structure 
of genèric Gaussian states of CV systems, and the origin of 
the promiscuity, are currently under way lf77ll . As an anticipa- 
tion, we can mention that promiscuity tends to survive even in 
the presence of mixedness of the state, but is destroyed by the 
loss of complete symmetry. The powerful consequences of 
the entanglement properties of GHZ/W states for experimen- 
tal implementations of CV quantum-information protocols are 
currently under investigation I78ll . 



VI. DECOHERENCE OF THREE-MODE STATES AND 
DECAY OF TRIPARTITE ENTANGLEMENT 

Remarkably, Gaussian states allow for a straightforward, 
analytical treatment of decoherence, accounting for the most 
common situations encountered in the experimental practice 
(like fibre propagations or cavity decays) and even for more 
general, 'exòtic' settings (like "squeezed" or common reser- 
voirs) 117911 . This agreeable feature, together with the possibil- 
ity - extensively exploited in this paper - of exactly computing 
several interesting benchmarks for such states, make Gaus- 
sian states a useful theoretical reference for investigating the 
effect of decoherence on the information and correlation con- 
tent of quantum states. Let us mention that the dissipative evo- 
lution of three-mode states has been considered in Ref. |80|. 
addressing SU(2,1) coherent states and focusing essentially 
on separability thresholds and telecloning efficiencies. In this 
section, we will explicitly show how the decoherence of three- 
mode Gaussian states may be exactly studied for any finite 
temperature, focusing on the evolution of the residual contan- 
gle as a measure of tripartite correlations. The results here 
obtained will be recovered in future work Il78ll . and applied 
to the study of the effect of decoherence on multiparty proto- 
cols of CV quantum communication with the classes of states 
we are addressing, thus completing the present analysis by in- 
vestigating its precise operational consequences. Concerning 
the general theory of open quantum dynamics, it is impossible 
here to give a detailed account of all the aspects of the Stan- 
dard theoretical frameworks. For an excel·lent critical review, 
focusing on the Standard treatment of open quantum systems 



in relation to quantum entanglement see Ref. 1 8 1 ] . In this am- 
ple review the authors discuss the importance of notions such 
as complete positivity, a physically motivated algebraic con- 
straint on the quantum dynamics, in relation to quantum en- 
tanglement, and analyze the entanglement power of heat baths 
versus their decohering properties. 

For continuous-variable systems, in the most customary 
and relevant instances the bath interacting with a set of n 
modes can be modeled by n independent continua of oscil- 
lators, coupled to the bath through a quadratic Hamiltonian 
Hi n t in rotating wave approximation, reading 



£ 

i=l 



v t (uj)[a\b t (uj) + a t b\(ui)\ du . 



(84) 



where bi(ui) stands for the annihilation operator of the ith 
continuum's mode labeled by the frequency w, whereas Ví(uj) 
represents the coupling of such a mode to the mode í of the 
system (assumed, for simplicity, to be real). The state of the 
bath is assumed to be stationary. Under the Born-Markov ap- 
proximation 118211 . the Hamiltonian Hi„ t leads, upon partial 
tracing over the bath, to the following master equation for the 
n modes of the system (in interaction picture) 18311 



^ 2 

i=l 



^,L[aí]g+(JV t + l)L[o,]g) 



(85) 



where the dot stands for time-derivative, the Lindblad super- 
operators are defined as L[ò]g = 2ògo^ — ó'òg — gò'ò, the 
couplings are 7.; = 2irvf(uji), whereas the coefficients Ni are 
defined in terms of the correlation functions (b\(üJi)bi(u)i)) = 
N^ where averages are computed over the state of the bath 
and lüí is the frequency of mode i. Notice that N is the num- 
ber of thermal photons present in the reservoir associated to 
mode i, related to the temperature Tj of the reservoir by the 
Bose statistics at null Chemical potential: 



Ni = 



kTi ■ 



1 



(86) 



In the derivation, we have also assumed (í>ï(wj)ò,(w,)) = 0, 
holding for a bath at thermal equilibrium. We will henceforth 
refer to a "homogeneous" bath in the case N = N and 7$ = 7 
for all í. 

Now, the master equation d85l > admits a simple and physi- 
cally transparent representation as a diffusion equation for the 
time-dependent characteristic function of the system í) 
H 



*«>*) = -£ 



(Xi Pi) 



dx q 
dpi 



) \Pi 



(87) 



where £ = (xi,pi, . . . ,x n ,p n ) is a phase-space vector, <j ioo = diag (2N + 1,2N^ + 1) and u> is the 2x2 sym- 
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plectic form [defined in Eq. Q]. The right hand side of the 
previous equation contains a deterministic drift term, which 
has the effect of damping the first moments to zero on a 
time scale of 7 /2 and a diffusion term with diffusion matrix 
(Tqo = ©™ =1 cr ioo . The essential point here is that Eq. d87i 
preserves the Gaussian character of the initial state, as can be 
straightforwardly checked for any initial CM cr by inserting 
the Gaussian characteristic function x(Ç, t) [see Eq. @] 

X(£,í) = e-if^^W^Í+^^í 

(where X are genèric initial first moments, <r(t) = I^ero + 
(1 - T^ctoo and T t = eie~ 7it / 2 l2) into the equation and 
verifying that it is indeed a solution. Notice that, for a homo- 
geneous bath, the diagonal matrices T t and (Jqo (providing a 
full characterisation of the bath) are both proportional to the 
identity. In order to keep track of the decay of correlations of 
Gaussian states, we are interested in the evolution of the ini- 
tial CM (To under the action of the bath which, recalling our 
previous Gaussian solution, is just described by 

o-(í) = It<r + (1 - r?)^ (88) 

This simple equation describes the dissipative evolution of the 
CM of any initial state under the action of a thermal envi- 
ronment and, at zero temperature, under the action of "pure 
losses" (recovered in the instance Ni — for i = 1, . . . , n). 
It yields a bàsic, significant example of 'Gaussian channeF, 
i.e. of a map mapping Gaussian states into Gaussian states un- 
der generally non unitary evolutions. Exploiting Eq. (18 8 1 and 
our previous findings, we can now study the exact evolution 
of the tripartite entanglement of Gaussian states under the de- 
coherent action of losses and thermal noise. For simplicity, 
we will mainly consider homogeneous baths. 

As a first general remark let us notice that, in the case of 
a zero temperature bath (N ~ 0), in which decoherence is 
entirely due to losses, the bipartite entanglement between any 
different partitions decays in time but persists for an infinite 
time. This is a general property of Gaussian entanglement 
ll79Íl under any many mode bipartition. The same fact is also 
true for the genuine tripartite entanglement, quantified by the 
residual contangle. If N 7^ 0, a finite time does exist for which 
tripartite quantum correlations disappear. In general, the two- 
mode entanglement between any given mode and any other of 
the remaining two modes vanishes before than the three-mode 
bipartite entanglement between such a mode and the other 
two [not surprisingly, as the former quantity is, at the begin- 
ning, bounded by the latter because of the CKW monogamy 
inequality (I47H. 

The main issue addressed in this analysis has consisted in 
inspecting the robustness of different forms of genuine tripar- 
tite entanglement, previously introduced in the paper. Notice 
that an analogous question has been addressed in the qubit 
scenario, by comparing the action of decoherence on the resid- 
ual tangle of the inequivalent sets of GHZ and W states: W 
states, which are by definition more robust under subsystem 
erasure, proved more robust under decoherence as well |84|. 
In our instance, the symmetric GHZ/W states constitute a 
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FIG. 5: (color online) Evolution of the residual Gaussian contangle 
G r T £S for GHZ/W states with local mixedness a = 2 (sòlid curves) 
and T states with local mixedness a = 2.8014 (dashed curves). Such 
states have equal initial residual contangle. The uppermost curves 
refer to a homogeneous bath with N = (pure losses), while the 
lowermost curves refer to a homogeneous bath with N = 1. As 
apparent, thermal photons are responsible for the vanishing of entan- 
glement at finite times. 



promising candidate for the role of most robust Gaussian tri- 
partite entangled states, as somehow expected. Evidence sup- 
porting this conjecture is shown in Fig.|5] where the evolution 
in different baths of the tripartite entanglement of GHZ/W 
states is compared to that of symmetric T states (at same ini- 
tial entanglement). No fully symmetric states with tripartite 
entanglement more robust than GHZ/W states were found by 
further numerical inspection. Qui te remarkably, the promiscu- 
ous sharing of quantum correlations, proper to GHZ/W states, 
appears to better preserve genuine multipartite entanglement 
against the action of decoherence. 

Notice also that, for a homogeneous bath and for all fully 
symmetric and bisymmetric three-mode states, the decoher- 
ence of the global bipartite entanglement of the state is 
the same as that of the corresponding equivalent two-mode 
states (obtained through unitary localization). Indeed, for any 
bisymmetric state which can be localized by an orthogonal 
transformation (like a beam-splitter), the unitary localization 
and the action of the decoherent map of Eq. J88i commute, 
because cr oo oc 1 is obviously preserved under orthogonal 
transformations (note that the bisymmetry of the state is main- 
tained through the channel, due to the symmetry of the latter). 
In such cases, the decoherence of the bipartite entanglement of 
the original three-mode state (with genuine tripartite correla- 
tions) is exactly equivalent to that of the corresponding initial 
two-mode state obtained by unitary localization. This equiva- 
lence breaks down, even for GHZ/W states which can be lo- 
calized through an (orthogonal) beam-splitter transformation, 
for non homogeneous baths, i.e. if the thermal photon num- 
bers Ni related to different modes are different [which is the 
case for different temperatures T, or for different freqüències 
üJí, according to Eq. J86H or if the couplings 7, are different. 
In this instance let us remark that unitary localization could 
provide a way to cope with decoherence, limiting its hinder- 
ing effect on entanglement. In fact, let us suppose that a given 
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amount of genuine tripartite entanglement is stored in a sym- 
metric (unitarily localizable) three-mode state and is meant to 
be exploited, at some (later) time, to implement tripartite pro- 
tocols. During the period going from its creation to its actual 
use such an entanglement decays under the action of deco- 
herence. Suppose the three modes involved in the process do 
not decay with the same rate (different 7;) or under the same 
amount of thermal photons (different Ni), then the obvious, 
optimal way to shield tripartite entanglement is concentrating 
it, by unitary localization, in the two least decoherent modes. 
The entanglement can then be redistributed among the three 
modes by a reversal unitary operation, just before employing 
the state. Of course, the concentration and distribution of en- 
tanglement require a high degree of non local control on two 
of the three-modes, which would not always be allowed in 
realistic operating conditions. 

The bipartite entanglement of GHZ/W states (under (1 + 
2)-mode bipartitions) decays slightly faster (in homogeneous 
baths with equal number of photons) than that of an initial 
pure two-mode squeezed state (also known as "twin-beam" 
state) with the same initial entanglement. In this respect, 
multimode entanglement is more fragile than two-mode, as 
the Hilbert space exposed to decoherence which contains it 
is larger. Notice that this claim does not refute the one of 
Ref. I80I1 . where SU(2,1) coherent states were found to be 
as robust as corresponding two-mode states, but only for the 
same total number of thermal photons in the multimode chan- 
nels. 



VII. CONCLUDING REMARKS AND OUTLOOK 

Gaussian states distinctively stand out in the infinite vari- 
ety of quantum states of continuous-variable systems, both for 
the analytic description they allow in terms of covariance ma- 
trices and symplectic operations, and for the high standards 
currently reached in their experimental production, manipu- 
lation and implementation for CV quantum information pro- 
cessing. Still, some recent results demonstrate that basically 
the current state of the art in the theoretical understanding and 
experimental control of CV entanglement is strongly push- 
ing towards the boundaries of the "ideal" realm of Gaussian 
states and Gaussian operations. For instance, Gaussian en- 
tanglement cannot be distilled by Gaussian operations alone 
I^l70ll7lll . and moreover Gaussian states are "extremal", in 
the sense that they are the least entangled among all states of 
CV systems with a given CM |85|. On the other hand, how- 
ever, some important pieces of knowledge in the theory of en- 
tanglement of Gaussian states are still lacking. The most im- 
portant asymptotic measures of entanglement endowed with 
a physical meaning, the entanglement cost and the entangle- 
ment of distillation cannot be computed, and the entanglement 
of formation is computable only in the special case of two- 
mode, symmetric Gaussian states l33ll . Moreover, when mov- 
ing to consider multipartite entanglement, many of the bàsic 
qüestions are still unanswered, much like in the case of multi- 
partite entanglement in states of many qubits. 

In this work we took a step ahead in the characterization 



of multipartite entanglement in Gaussian states. We focused 
on the prototypical structure of a CV system with more than 
two parties, that is a three-mode system prepared in a Gaus- 
sian state. We completed the elegant qualificative classifica- 
tion of separability in three-mode Gaussian states provided 
in Ref. fïlll with an exhaustive, quantitative characterization 
of the various forms of quantum correlations that can arise 
among the three parties. We then exploited some recent re- 
sults on entanglement sharing in multimode Gaussian states 
fÏ2ll that prové that CV entanglement in these states is indeed 
monogamous in the sense of the Coffman-Kundu-Wootters 
monogamy inequality 0|. We next defined a measure of gen- 
uine tripartite entanglement, the residual continuous-variable 
tangle, that turns out to be an entanglement monotone under 
tripartite Gaussian LOCC lfÏ2ll . 

We started our analysis by giving a complete characteriza- 
tion of pure and mixed three-mode Gaussian states, and de- 
riving the Standard forms of the covariance matrices that are 
similar to those known for two-mode states 14111 . In particu- 
lar, a genèric pure three-mode Gaussian states is completely 
specified, in Standard form, by three parameters, which are the 
purities (determinants of the CMs) of the reduced states for 
each mode. We determined analytically the general expres- 
sion of the genuine tripartite entanglement in pure three-mode 
Gaussian states, and studied its properties in comparison with 
the bipartite entanglement across different partitions. We in- 
vestigated the sharing structure underlying the distribution of 
quantum correlations among three modes in arbitrary Gaus- 
sian states, much on the same lines as those followed in the 
case of states of three qubits Ió4ll . 

Remarkably, we found a completely unique feature, namely 
that that there exists a special class of states, the pure, symmet- 
ric, three-mode squeezed states, which simultaneously maxi- 
mize the genuine tripartite entanglement and the bipartite en- 
tanglement in the reduced states of any pair of modes. This 
property, which has no counterpart in finite-dimensional sys- 
tems, can be understood as the promiscuous sharing of CV en- 
tanglement. The states exhibiting this peculiar sharing struc- 
ture, named CV "GHZ/W" states for self-explaining reasons, 
are automatically preselected as optimal carriers of quantum 
information over lossy channels, and we proved that they in- 
deed are. In fact, we concluded our work with a detailed anal- 
ysis of the effects of decoherence on three-mode Gaussian 
states and the decay of tripartite entanglement. This study 
yielded that the GHZ/W states are the most robust three-party 
entangled Gaussian states against decoherence. 

We believe that the collection of results presented here, al- 
though remarkable on its own, is however only the tip of an 
iceberg. Three-mode Gaussian states, the perfect test-ground 
for the understanding of some genèric traits of multipartite 
entanglement in CV systems, need to be analyzed in a deeper 
future perspective. This primarily includes the characteriza- 
tion of those classes of tripartite entangled states with pecu- 
liar properties, with a particular care towards their state en- 
gineering in quantum optical settings. This analysis is cur- 
rently under way 17711 . The (closely related) usefulness of such 
states for existing and maybe novel protocols of CV quantum 
communication, able to take advantage from the promiscuous 
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sharing, is also being investigated 1 78 ] . 

From a broader theoretical standpoint, further research 
stemming from the present work should probably be directed 
along two main directions. The first one concerns proving a 
general monogamy inequality in all multimode states of CV 
systems, in analogy to what has been recently established for 
arbitrary states of multiqubit systems |4]. Such a proof would 
then lead to a multimode generalization of the residual con- 
tangle. The second, long-term direction is the investigation 
of the qualitative and quantitative aspects of entanglement in 
genèric non-Gaussian states of CV systems. In this context, 
singling out exòtic states with enhanced promiscuous sharing 
of quantum correlations and with a monogamy of entangle- 



ment stretched to its limits, appears as an exciting perspec- 
tive, and might open very promising perspectives for the ma- 
nipulation, transfer, and control of quantum information with 
continuous variables. 
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